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Abstract 

We show the equivalence between the three approximation schemes for self- 
interacting (1+1)-D scalar field theories. Based on rigorous results of [1, 2], we 
are able to prove that the Gaussian approximation is very precise for certain 
limits of coupling constants. The A0 4 + acp 2 model will be used as a concrete 
application. 

I. INTRODUCTION 

We start by clarifying the terms which appear in the title. For a given theory of self- 
interacting quantum fields, by Gaussian approximation we mean the Gaussian part of the 
interacting measure. By mean field approximation we understand the leading term of the 
expansion given in [1, 2] and by variational approximation we understand the old variational 
technique performed with some particular trial states [3-5]. 

The novelty of this paper is the method of extracting the Gaussian peace of the interacting 
measure. Also we hope that the equivalence between the three methods will lead to a better 
understanding of the self-interacting field theories. The expansion around the mean field 
(MF) approximation have been successfully used in [1, 2] to A -2 Y2n=i a n i^) 2 ™ interactions, 
for small A. The key role of this approximation is that, after the translation to — <Pmf, it 
brings all the coupling constants near to zero where the cluster expansion [6] can be applied. 



The results show that the mean field approximation incorporates the nonanalytic part of the 
Schwinger functions. Our Gaussian approximation will play the same role. For A — > 0, oo, 
it provides a transformation of the field which brings all the coupling constants near to zero. 
Actually, for A — > 0, the Gaussian approximation reduces to the mean field approximation. 
A similar procedure have been used by Glimm et al in [7] to prove the equivalence between 
the model :A0 4 + <t0 2 :, with large A, and the model :A0 4 + <70 2 : with small A and negative a. 
However, the variational approximation scheme (especially in the form of [5]) became very 
popular among physicists even thought that there were no estimates of the errors related to 
this scheme. The equivalence of the three methods allows us to use the exact results of [1, 
2] to estimate these errors. 



II. OUTLINE OF THE STRATEGY 

We consider in this paper (1+1)-D scalar fields with self-interactions. To build such theories, 
one can follow the strategy presented in [8] . Consider first a space cut-off interaction: 

U(s)= ( rf 2 x:\/(0(x)): mo , (1) 

J\x\<s 

where the normal ordering is with respect to the vacuum of the free field of mass m . The 
next step is the investigation of the cut-off interacting measure: 

where /w is the Gaussian measure corresponding to the covariance C mn — — : tt and 

— A + 771q 

mean £. Note that the additional term to the potential U cancels the mass and the mean 
of the measure // ,£- The last step is the investigation of the (thermodynamic) limit of 
this measure: lim /i s . In the case when the limit is well defined, one can build the physical 
Minkowsky field by using one of the well known reconstruction methods. The hardest part of 
such a program is obvious the thermodynamic limit. Our strategy of pursuing this program 
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will be as follow. The investigation of fi s is equivalent with the investigation of the moments 
(Schwinger functions): 

(<j){xx) ...<j){x n )) s = j (j){x 1 ) ...(j>(x n )dii s . (3) 
The thermodynamic limit can be studied by considering the dynamical system: 



d 
ds 



1 * X 



X, (4) 



(cj)(x 1 )(j)(x 2 )) s 

v - I 

where X is an infinite vector of whom expression can be deduced by simply taking the 
derivative of the right part of ||] in respect to s. The thermodynamic limit is reduced now to 
the study of the stable fixed points of the dynamical system |4|: if the initial conditions for 
this dynamical system lie in the basin of attraction of a stable fixed point then we found the 
thermodynamic limit. We will show that the initial conditions can be modified by changing 
the boundary conditions of the field. Actually, in the Gaussian approximation we can choose 
the initial conditions to superimpose over the stable fixed points. This will lead to a set of 
self-consistency equations which will provide us an approximative value for the expectation 
value and the mass of the field. For interactions of the type: A~ 2 P (A0), the solution of this 
set of equations reveals two things: in the limit A — ► we recover what Glimm et al call 
the mean field approximation. In the limit A — > oo (strong couplings), the mass increases 
much faster than A so by a rescaling of the field we are again in the small coupling constants 
domain. To show the equivalence between Gaussian and Variational approximations, we 
compute the vacuum energy in Gaussian approximation and show that it coincides with the 
expression given in [5]. 



III. THE GAUSSIAN APPROXIMATION 



Let us start by specifying our notations. In general, by : • : M we denote the normal ordering 
with respect to some measure \x. If the measure corresponds to a covariance C and mean £, 
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we denote it by : • :c,£ and further, if the covariance corresponds to some mass, m, then we 
use : ■ : m £. Instead of : • we will use : • : s . We denote the expectation values with respect 
to some measure fi by (■) and we apply the same shorthands as for the normal ordering. 

A. The dynamics of the expectation values 

Starting from: 

... , , ( u JHpl) -4(x n ) e- u 'Wd(Mo 

(0 (-1) -4 (*„)>. = /e -^ %Q > (5) 



then: 



^ s (cj ) {x 1 )..4{xn)) s =^j^_ dx [{V'(<j>(x))) s -V'(<f>(x))} ( f>(x 1 )... ( j>(xn)^ (6) 
or, shortly: 

j- g (0 ..4 (x n )) s = ([([/' (8 r )) s - U> (Sr)] <f> (x,) ..4 (x n )) s , (7) 

where T is the curve \x\ — s and 5 is the Dirac delta function. As is stated now, the above 
equation has the initial condition: 

((p (xi) ..4 (x n )) s=0 = (0 (xt) ..4 (x n )) moif , (8) 

in particular: (<p (x)) s=0 = £ and (: 4> (x) : s=0 : <p (y) : s=0 ) = C mo (x - y). 

B. The change of covariance 

We can allow us a little more liberty in choosing the initial conditions. If one uses the 
Gaussian perturbation identity 0: 

d/J, (c - 1+v) -i = Z' 1 exp (- : V : c ) d/J, c , (9) 

where : V \c= - J v (x) : <p (x) 2 and Z = J exp (— : V :c) d<f>c, it follows that, up to a 



change of boundary conditions: 
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where: 



['(,)= / d 2 xV(<P(x))= I d 2 x {:V(<j ) (x)): mo -^:(<f ) (x)-0 2 -, 

'\x\<s ^ |x|<s 



:n) 



and dun f is the Gaussian measure corresponding to the covariance C m = — and 

—A + m z 

mean £. Now it is the time to discuss the boundary conditions. One can see that from very 
beginning some boundary conditions were imposed. For example in Q the additional term 
cancels the mass and mean of the measure dfio,£ only for \x\ < s. This cut-off is equivalent 
with adding an interaction in the region \x\ > s. In the limit s — ► oo, it provides us the 
external conditions. One can read [9] for a discussion of this boundary conditions. The same 
remarks for the Gaussian perturbation identity was cut off, the effect being the change 
of the external conditions. In the light of this remarks, one can see that all the time when m 
and £ are changed, the boundary conditions are changed or, reciprocally, we can change m 
and £ by a modification of the boundary conditions. As a result, we can control the initial 
conditions of the dynamical system ^ by imposing different boundary conditions: 

- (0 ( Xl ) ...0 (*„)), = ( [(U (5 T )) s - U (5 T )\ (x x ) ...0 (x n )) ^ 

(12) 

(0 (Xl) ...0 (X n )) s=0 = (0 (Xl) ...(f) (X n )) m$ 

If the system is multiphasic, then O will have more than one stable fixed point. The above 
procedure can be used in the following way to select a pure phase: enforce the thermody- 
namic limit to converge to a fixed stable point by choosing the initial conditions in the basin 
of attraction of this point. There is now a nice picture of critical phenomena in terms of the 
bifurcation theory. Suppose [12] has a degenerate fixed point and there is a bifurcation point 
where the degeneracy is lifted. Then the critical values of the parameters are given by the 
bifurcation point. 
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C. The Gaussian approximation 



We consider that fi s is approximately Gaussian and symmetric at translations, or one can 
think that from now we are interested only in the symmetric Gaussian part (and the mean) 
of the measure fi s . The major property of a Gaussian measure is that the normal ordered 
powers are orthogonal: 

(: r : s : m : s ) s = if m ^ n. (13) 

We will select the following two equations of the dynamical system which in the Gaussian 
approximation write: 

' ^ (0 (0)) s = -T M (V (0 (x)j) (: <f> (6 r ) :.: <f> (0) : s ) s 

< 

k ^ (<: 0(0) :?} s - (: 4>(*) 4o) mo ) = (v{4>W)) < : ^ (fr) 0(0) (14) 

where T UiS (■) means the coefficient of : <p n : s of the given expression. The point x satisfies 
\x\ = s and we have used the symmetry of the problem to separate the terms. Note that the 
second equation can be written in the above form because (tj (Sr)j —U (Sr) has no free term 
and : (ft 2 : s — : (ft '? 8 is a constant. We have taken the combination (: (p (0) 's) a — \' (ft (x) : m ) mo 
because it is a finite quantity and for the reason explained in the following. After one changes 
the normal ordering in V (0) from : ■ : mo to : • : s , 

:0(*r :mo=limELo^X 

a — >0 

x £i { exp (y ( (: (x) ~ < : (x) i>0 + a (0 } : Hx)k - ' (15) 

relation which is valid in the Gaussian approximation, one can see that T n s (v {<p{x))^ is 
a function only of the combination (: </>(x) :^) s — (: <f>(x) : m ) mo an< ^ the expectation value 
(4> {x)) s . Now, because we consider also the symmetric part of fi s at translations, the above 
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system can be written in the form: 



d 



- (0 (x)) s = -T M ( V (0 (*)) ) (: (S r ) :.: (0) :« 



"to 



-T 2 , s (v (*(*))) (: 2 (^r) : S :0(O) :^ 



(16) 



Let us rewrite this system in the following notations: 



(17) 



m 



Then: 



( d ~§ = -Ti,s (X,Y) (:0(<5 r ) : S :0(O): 



< ^ = -T 2iS (X,F)(:0 2 (5 r ): s :0(O) 



: 2 ) 

s/ s 



(18) 



1 



X(0)={,K(0) = -h.(^}. 

It is easily now to find the fixed points of the dynamical system (more exactly, the first two 
coordinates of the fixed points), which are given by the conditions: 



T 1>3 (X,Y) = 0,T 2>9 (X,Y)=0. 



(19) 



Because we do not have to much information about the basin of attraction of the fixed 
points, the best thing we can do is to superimpose the initial conditions over the stable 
points. This leads to our self-consistency equations: 



rrtc 



Tl ' s ^'4^ ln (^ ))=0 ' T2 ' s(e '- ln ' V 



1 , / m\ 
An \ m 2 



0, 



(20) 



which provide us the Gaussian approximation of the expectation value (0 (x))^^ = £ and 
the two point Schwinger function (: 0(zi) : s ~>oo- <j>{x 2 ) : S - +O o) s ^ o = C m ( x i, x 2)- 



IV. THE EQUIVALENCE WITH THE VARIATIONAL METHOD 



We calculate the energy of the vacuum state in the Gaussian approximation. Will follow 
that our expression is the same with that of [5]. Moreover, our self-consistency equations 
can be derived by minimizing this energy with respect to £ and m. 

The term Z, which appears in the Gaussian perturbation identity, was unimportant for the 
expectation values but will be essential for the calculus of the vacuum energy According to 
||, the energy per unit length can be calculated as: 

lnje- u '^d m 



e = — lim 

s— >oo 



TCS A 



(21) 



After the change of the covariance: 



lim 

s— >oo 



InZ-ife-^dfic^ 



7TS Z 



(22) 



with: 



m — m, 



cxp 



which can be calculated formally || as: 



1, 



d 2 x : (0 (x) - 



■m 



dfi, 



o,5- 



(23) 



Z = exp -Tr < In 1 + 

1 2 I V -A + ml 



m 



m — m. 



(24) 



-A + ml 

In this form, the above quantity is infinite. Taking into account that when we have performed 



the change of covariance we have considered the integral J d' 



x 



\x\ < s, we should do the same thing for Z. Considering the operator: 



K = In 1 



m 



-A + ml 



m 



- A + ml ' 



{x) ~ -mo only on 



(25) 



with the kernel: 

K(x,y) 

it follows: 



d 2 k { In 1 + 



1 m 2 



1 m 



47r 2 kh2 + mlJ An 2 k 2 + ml 



-ik(x-y) 



(26) 



Tr 



K 



x\<s 



d 2 x ft In 1 



1 2 2 

1 m — 777,0 

47r 2 kh2 + ml 



1 2 2 

1 777 — 

47T 2 k 2 + ml 



(27) 



if we restrict the domain of integration at \x\ < s. Finally: 



Tr 



K 



US 2 ( 2 2 2 m 2 

m — m n — mm — ^ 

4tt V ml 



(28) 



In the Gaussian approximation, the other term of e: - — — is given by T 0jS ( V ) . 



7TS Z 



In consequence: 



/ ~ \ 1 / 2 2 2 i m ^ 

e = T , s K + — m -m -m In — = 



(29) 



The last expression is completely equivalent with that of [5]. Indeed, for an analytic function, 
F, it follows from 1131 that: 



To, [:F(<f>(x)) : mo ] = E~ 



1 <9 n F(a;) 



n! <9x r 



x=0 



<9 n 

<9a r ' 



a , m 



exp — In — - + a£ 



57T m z 



a=0 



n=0 2^fc=0 



(n — k)\k\ dx n 



tn—k 



x=0 



^7r m 2 



<9 fc exp 



<9a fc 



a=0 



Eoo \-~\n 
n=0 2^k=0 



d n F 



(n — k)\k\ dx r 



tn—k 



x=0 



Qk 

da k 



a=0 



1 d n F 
l^n=o 2^=o ( n _ 2p )!(2p)! 



x=0 



p! \87r m 2 



2\ P 



1 a n F 



2p)! <9x n 



jn— 2p 



x=0 



exp (±ln™° — ) F 

\8tt m 2 dx 2 ) 



x 



With this result we have: 



e (m, £) = exp I — In 



1 , ml d 2 



In m 2 dx 2 



m 



v{t)- T {x-ty 



I / ™2 

H ( m 2 — m„ — m 2 In 

8vr V ™ 



} - exp ( i_ l n ^ V (0 + i K - O ■ (30) 



87r m 2 d£ 



57T 



We show in the following that our self-consistency equations can be also derived by mini- 
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mizing e in respect to £ and m. Indeed: 

1 d n V(x) 



: V (4>(x)) : mo =^^=0 



n\ dx n 



x=0 



a=0 



Eoo 
n- 



i <w 



n! <9x n 



x=0 



<9 n f / a 2 m 2 . 



a=0 



(31) 



d n f /a 2 , m 2 A] 



a=0 = ae" 



x=0 

2 

Using the same notation as before, Y = gMn^jf, for the second equations we have: 



T, 



2,s 



■ V {<P (X)) -mo =E~=0 



1 d n V 



n\ dx r - 



n (n — 1) 







,n-2 



da n ~ 2 



x=0 



. a , mZ 
exp — in — ^ + at, 

7r rrr 



a=0 



d n V 


<9 n 


dx n 






x=0 




( d n V 







{a 2 exp (a 2 F + a£)} 



a=0 



(9™ 

9a 7 



{exp (a 2 F + <)} 



x=0 



a=0, 



9 / d n V 



dY \ dx r 



i=0/ 



9a™ 



{exp (a 2 F + a£)} 



a=0 



(32) 



_d_ 

dY 



To, s ( V 



2 1 9 / d H 



z "" =0 n! <9F 



^ (x - S 



dx n 



Using that: 

_d_ ■ 
dY . 



Tn, V 



w{ £ - m2 (i + Y 



Qn 

da n 



{exp (a 2 Y + of)} 



x=0/ 



a=0 



and expanding the second term of |32], it follows: 

":y(0(x)): mo 



■2,8 



de_ 
dY' 



(33) 



(34) 
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We have proven that the Gaussian approximation of the vacuum energy density is equal with 
the expression obtained in [5] and that our self-consistency equations can be also derived by 
minimizing the energy density with respect to £ and m. In consequence the two methods 
are completely equivalent. 



V. THE V 



\<j) 4 + a<f> 2 MODEL 



We apply in this section the Gaussian approximation to this particular model. The self- 
consistency equations are: 



' T ltS \:V(4>) : m J =i(3Aln^| + 47rA£ 2 + 27r<7) =0 



7T 



m 



(35) 



2,s 



3A / rrin A m 2 



■V{4>): mo =— ln^ + 47r£ 2 — = 0. 

7r V m z I 2 



Note that the first equation always has the solution £ = 0. Let us analyze first this situation. 
The second equation leads to: 



m 



3 A / mhx ( 2ncr 
W ' 1 



7T 



3A 



exp 



V 3A 



(36) 



where W is the Lambert W function of rank zero. For a > 0, we can introduce the classical 
mass: ml = V" (0) = 2a. One can see that, if we choose m = m c , then m c is a solution of 
the above equation. For a < 0, the solution £ = becomes unstable. The other solutions 
must satisfy: 

'2tt 



m 2 = 8A£ 2 . 



(37) 



For a < we can define the mean field approximation: (£ c ) 2 = — — and m 2 = V" (£ c ) 

2A 



8AC 2 



-4(7. One can see that, if we start with m = m c , then the mean field approximation 



is a solution of our self-consistency equations. This is true for more general interactions and 
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it shows the link between the Gaussian and mean field approximations. Fortunately this is 
not the only solution. The general solution of |37] (for tjiq = m c ) is: 



^2 — 3 TIr / 2lIO / 2-ko 



where can be the Lambert function of rank or —1. It follows that: 

2na 



(38) 



' -3 /2na ( 2na 
W-i \ "ttv exp 



, c2 



47T 



3A 



3A 



< 



3 /27TCT /2tTCT 

^hr exp br 



3A 

2vra 
~3A 



> -1 



(39) 



is the solution which coincides with £ c . It gives the leading term of ((f)) in the limit A — ► 0. 
The other solution is: 

2txo 



— 3 TTr / 2txo ( 2na , , 



-3 Trr /27ra /27ra\\ 2vra 



(40) 



which gives the leading term of ((f)) in the limit A — > oo. Moreover, even when the classical 
picture is lost, a > 0, in which case the potential does show only one minima, there is a 
solution which shows symmetry breaking in the limit A — > oo: 



e 2 



-3 

47T 



6A 



exp 



27TCT 

"3A~ 



(41) 



The nice behavior of these solutions relies on the fact that, for A — > 0, the first solution 
keep the mass constant (while the coupling constants goes to zero) and, for A — ► oo, the 
second solution gives a very large mass such that the ratio between the coupling constants 
and mass goes to zero. To see this, let us calculate the expression of the potential after we 
apply the Gaussian transformations. It follows from the analysis of the last sections that 
the interacting measure is equal to (up to our boundary conditions): 



C.-ii 



(42) 
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and after the change of the normal ordering in the exponent, with the values of m and £ 
given by the self-consistency equations, the coefficients of :0 2 : and :</>: cancel out, the result 
being: 

dfjiA = VW**(*) 4 ^*(*)Ve^ Cmjf . (43) 
Finally we rescale the mass at the unity by using the rescaling identity ||: 

dfiAf = Zl}e- S ^' l W^ + k^' ) ^diJL (hAt (44) 



with A' = A/m 2 . Note that we have omitted the free term of the potential, e (m, £), which 
only shifts the energies and is unimportant for the expectations values. Thus one can see 
that the small coupling constant regime is achieved for large values of £. However, even 
for large values of £, the convergence of the thermodynamic limit cannot be proved by an 
ordinary cluster expansion. This is due to the fact that the potential V = g|2</> 4 + <^0 3 is not 
bounded from below uniformly in £. When the field is localized near to — £, it behaves as 
— £ 2 . Nevertheless, an expansion in phase boundaries [1] will solve the problem. We can give 
now the asymptotic expansion of Schwinger functions. As usual, we will apply the formula 



for integration by parts to the interaction 03 



f:0(x) m :R (0) d/i A = JJ Ci (x - y) :<f> (x) m 1 : j ^|^y - R J^^ } dydfx A . 

(45) 

The formula is still true for the case when the normal ordering is with respect to a Gaussian 
measure with a nonzero mean value (see Appendix). For (0) we have to apply this formula 
three times to find the first correction to £. Redenoting A' as A, it follows: 

((f) (x)) = J (x) dfi A = /(:</> (z) +0^A 

(46) 

= Z+wUC 1 (x) 3 dx-lfd (x) C x (y) C, (x - yf dxdy] + o (l/£ 4 ) . 
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For the connected part of the two point Schwinger function, it follows: 



(0 (*) (y)> - (0 (*)> (0 (y)> = / : (x) : u : (y) : U <fc A + o (l/£ 3 ) 



(47) 



Ci (x-y) + ^ r //Ci(a;-z)C7i( 2 /-«) d u)^ dzdu -\- o (l/£^) . 



If we go back to the original scale, we can conclude: 



<</>(*)>=£ + 



0.021 



< 



(48) 



(0 (:r) (y)) T = C m (x - y) + 



5.6 x 10" 4 



+o(i/e) 



VI. CONCLUSIONS 



There are two important consequences of the analysis. First, it was unclear before what 
is the meaning of the parameter m. Second, it follows from the last section that this 
approximation is very precise for certain range of coupling constants. In that regime, it is 
easyly to see that m is in fact an approximation of the self-interacting field mass, i.e. the 



see that it is completely determined by only one parameter, the mean value of the field, 
which is experimentally measurable. 
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About the A0 4 + <r0 2 model, one can 
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VII. APPENDIX 



We prove first that the Wick powers are orthogonally if and only if the measure is Gaussian 



i.e.: {e a<t> ) = exp ^ (: : 2 ) + a (0) 



. This can be seen from: 



(e^) = (0e a< ^) = <: : e a<t> ) + (0) <0e a< ^) = <: :: e a< ^ :) (e a<t> ) + (0) <e a ^) 



Using the orthogonality of the Wick powers we can continue 



d 



(49) 



-<e^)=a<:0: 2 ><e^) + <0> 



(50) 



equation which integrated out (together with the initial condition (e a *)| a=0 = 1) leads to 
the desired expression. The other implication follows from 



exp(a0 + &v?) ( , . , v (exp(a0 + &y?)) 

: exp (a0) :: exp (6(/9j := - — - — — — - - v =: exp (a<p + hp) : 



(exp a0) (exp 699) 



(exp a0) (exp 69?) 



(51) 



so, for and </? jointly Gaussian variables (with mean different by zero): 
(: exp (a0) :: exp {bp) :) 

exp (| (: a0 + 6^ : 2 ) + (a0 + fry?)) 



(52) 



exp 



^(:0: 2 > + a(0> 



exp 



§(:y?: 2 >+Mv?> 



j = exp (a& (: (f) :: ip :)) 



Taking derivatives in respect to a and b and then the limit a, b — > 0, we can form any Wick 
power of and </?. Because of the expression on the right side, it is obvious that the Wick 
powers are orthogonally. 
Integration by parts (first formulation) 



/ : (/) -c^ R (</>) duct = J f(x)C(x-y) 8R/8<f> (y) dxdyd\x CA 



(53) 
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It is enough to prove this formula for the case: R (0) = exp [i<j> (g)]. In this case: 
/ : {f) e^dnc* = / : e<W : ^ 



. d f , . d exp [-JC ( 5 + \f,g + \f) + i(^g + A/)] 



'dA J ( e W)> 



dA exp[-iC(A/,A/)+i(e,A/)] 



iC (/, 5) exp [-±C (<?, (?) + i (&</)]= i / / (*) C(x-y)g (y) dxdy (e^) ( 54 ) 



r t( \r<( \ 6 exp ^ 
= J f(x)C{x-y) dxdydncz, 

where (•, •) is the scalar product in L 2 (M 2 ) and / and g were considered in L x (M 2 ). Next 
we show that: 5 : (f) (/)" : /5<f)(x) =: 5<f)(f) n /5<f)(x) :, where the normal ordering is with 
respect to an arbitrary measure, which is helpful when one computes SR/S<f>. Indeed: 



5 5 d n 

— :o{f)": -—^-:,^f'- 



5(f) (x) 



5(f) (x) d\ r ' 



5 d n e x ^ 



=0 5(f) (x) dA B (eW) 



A-0 



d n Xf(x)e 



d\n ( e A0(/)) 



n/ (x) 



d n-l e XHf) 



A=0 



dA n-l ( e A<H/)) 



n/ (s) : (/) 



71-1 



A=0 



(55) 



;0(/) n :. 



' 5(f) (x) 

Integration by parts (second formulation) 



/ : (/)" '-c,a R (0) due* = I f ( x ) c ( x - y) ■ UT -u $ R / S( f> (y) dxdyd/i c ^ 



(56) 



We start with the identity: 



: i./T :c*=: (/) 'c* (/)"" 1 - f(x)C(x-y): 



. 50 (/) 



71-1 



50 (y) 



dxdy, 



(57) 
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which can be proven as follow: 



HIT :<*=^:eW) 



d n exp[A0(/)] 



A=0 



rfA" (exp [A0(/)]> 



A=0 



exp [\<f> (/)- |C7 (/,/)- A (£,/)] 



A=0 



dA 



n-l 



[*(/)- (6 /)] t^M -AC(/,/) 6Xp[A0(/)] 



(exp [A0(/)]> 



Jn— 1 



=: (/) : r . t : (/)"" -_//(*) C (s - y) 



(exp [A0(/)]>^ 
5 exp[A0(/)] 



A=0 



(58) 



50 (y) (exp [A0(/)]>^ 



dxdy 



A=0 



=: (/) :c&: (/)" ' -Jf(x)C(x-y) — — : (/)" ' : c>€ dsdy. 



50 (y) 



Then we can continue: 



/ : !/)" :c, f (0) d»cA = I ■ (/) ^: (/)" ! ^ A (0) d/** 



JSSf(x)C(x-y) : (Z)"" 1 :ca R (0) ^dyd/^ 



= fff(x)C(x-y) [: : c>€ i2 (0)] cfadyd^ 



(59) 



1 1 1 f( x ) c ( x -y) 7-7-T : (/)" # (0) dxdydnc^ 



50 (y) 

which leads to the desired expression. 



17 



REFERENCES 

[1] Glimm J, Jaffe A, Spencer T: Ann. Phys. 101, 610 (1976) 
[2] Imbrie J Z, Commun. Math. Phys. 82, 261 (1981) 

[3] Stevenson P M, Phys. Rev. D 30, 1712 (1984); ibid 32, 1389 (1985); ibid 33, 2305 (1986); 
[4] Ingermanson R, Nucl. Phys. B266, 620 (1986) 

[5] Coleman S, Aspects of Symmetry, Ch. 6, Cambridge Univ. Press, (1994-1996) 

[6] Glimm J, Jaffe A, Quantum Physics, Berlin, Heilderberg, New York: Springer (1981) 

[7] Glimm J, Jaffe A, Spencer T: Existence of Phase Transitions for QF, Marseille Conference 
(1975) 

[8] Simon B: The P(<fi) 2 Euclidean (Quantum) Field Theory, Princeton: Princeton Univer- 
sity Press (1974) 

[9] Frohlich J, in Invariant Wave Equation, Proc. of Eltore Majorana (1978) 



18 



